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Abstract 

Passing current at given threshold voltages through a metal/insulator/metal sandwich structure 
device may change its resistive state. Such resistive switching is unique to nanoscale devices, but its 
underlying physical mechanism remains unknown. We show that the different resistive states are 
due to different spontaneously charged states, characterized by different ‘band bending’ solutions 
of Poisson’s equation. For an insulator with mainly donor type defects, the low-resistivity state is 
characterized by a negatively charged insulator due to convex band bending, and the high-resistivity 
state by a positively charged insulator due to concave band bending; vice versa for insulators with 
mainly acceptor type defects. These multiple solutions coexist only for nanoscale devices and for 
bias voltages limited by the switching threshold values, where the system charge spontaneously 
changes and the system switches to another resistive state. We outline the general principles how 
this functionality depends on material properties and defect abundance of the insulator ‘storage 
medium’, and propose a new magnetic memristor device with increased storage capacity. 
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I. INTRODUCTION 


The existence of resistive switching memories, also known as memristors or Resistive 
Random Access Memories (ReRAM or RRAM), was predicted from circuit theory nearly 
half a century ago [l, 2| and functional devices have been demonstrated within the past 
decade [3l-l7|. The prototypical memristor is composed of a thin slab of Ti02 sandwiched 
between two metal electrodes js], but memristors have been realized using a wide range of 
insulator materials as storage medium jb, 7|. By applying bias voltages in excess of given 
threshold values, the resistivity of the storage medium can be changed from a low-resistivity 
state (LRS) to a high-resistivity state (HRS) and back, over a large number of cycles j^. A 


revolution in information technology is only hinder^ 
physical principles that underlie resistive switching 


by the lack of understanding of the 


mla, 


lol. 


Thus far, two types of mechanisms have been suggested for this resistive switching. Elec¬ 
tronic mechanisms describe the switching as electrons being trapped and un-trapped inside 
the insulator storage medium, such that the differently charged systems would exhibit dif¬ 
ferent resistivities UHlSj]. Ion drift mechanisms, on the other hand, attribute the switching 
to a concerted diffusion of atoms/ions that leads to a controlled and reproducible 


brmation 


and breaking of conductive hlaments through the insulator storage medium 14j-l22l|. The 
proposed electronic mechanisms suffer from the voltage-time dilemma jo] because the average 
lifetimes for the asserted metastable electron-trapping states are much shorter than typical 
memristor operation times, and ensueingly, such a system would spontaneously forget which 
resistive state it is in. Ion drift mechanisms offer arbitrarily long lifetimes for the different 
resistivity states, because the corresponding hlament structures (formed or broken) are as¬ 
sumed to be stable at room temperature. This stability, however, implies that the formation 



predict switching times in the order of 10 ^ s. Moreover, neither mechanism justifies typi¬ 
cal on/off resistivity ratios (> 10^ for Ti 02 a. a„d even .neh la^ge. fo. othe. n,ate.a.s). 


For example, assuming the ion drift mechanism for Ti02, conduc 


Magneli phase hlaments increases only by a factor of 10^-10^ 8l. 1201 . 


Ivity through observed 


We present a new electronic mechanism for resistive switching. The LRS, HRS, and 
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possible intermediate resistivity states (IRS) correspond to different ‘band bending’ solntions 
of Poisson’s equation, which describes the electrostatic potential and charge stored inside 
the insulator storage medium. This implies that by controlling the defect abundance and 
size of the insulator material, one can construct a resistive memory using arbitrary insulator 
materials sandwiched between suitable metal electrodes. The present theory of resistive 
switching is universal, and there is no need to invoke different types of switching mechanisms 


for resistive memories made from different types of materials 
specihc correlated materials 


28 


or specihc theories for e.g. 


29| . The variationally calculated electrostatics, or band 


bending solutions presented here are general for any metal/insulator/metal structures (not 
only memristors), and as we show, in stark contrast to the typical approximations made in 
previous theories of resistive switching. The present theory does not exclude the existence 

3J3, 


of the formation and breaking of conductive hlaments during memristor operation 


but we propose that it has a different role in resistive switching, as was assumed previously. 

This work is organized as follows. First, the electrostatic model for different resistive 
states is presented, followed by a model device and an electrostatic theory for resistive 
switching. We then discuss the inhuence of material parameters on the multiple resistive 
states, and outline the design principles how to craft memristors using arbitrary insulator 
materials as storage medium, followed by a discussion of the role of storage media hlamen- 
tation in resistive switching. Finally, we propose a new magnetic memristor device with 
increased storage capacity, and summarize our results in the Conclusion. Details of our 
variational calculations are given in the Appendix. 


II. ELECTROSTATIC DESCRIPTION 


The key point that is missing in previous switching mechanism models—electronic or ion 
drift based—is that they all have assumed the electrostatic potential (j){z) within the insulator 
storage medium to exhibit either typical metal-insulator band bendings, or no band bending 
at all {q, ^ (here e0 marks the position of the conduction band minimum, see App. lAl) . 
Such band bendings follow approximate (analytic) solutions to the one-dimensional Poisson’s 
equation 


d‘^(j) p{z) 

dz"^ eeo 


( 1 ) 
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along z perpendicular to the metal-insulator interface, where e and eo are the relative and 
vacuum permittivities, respectively, and the charge density p{z) is approximated by a step 
function. This approximation for p{z) is justihed for a semiconductor with a small concen¬ 
tration of (shallow) dopants, but not necessarily in memristors, where defect concentrations 
can be 10^^ cm“^ (= 1 nm“^) or even larger |8|. In fact, p is a functional that depends on 
both 2 ; and (l>{z), i.e., p = p[ 0 (z), 2 ;], and must be solved self-consistently together with (j){z). 
We thus solve Eq. ([1]) variationally by minimizing the functional 

FMz)] = £ dz + 4'£) ■ (2) 

Here a and b denote the z coordinates of the two insulator/metal interfaces, and 0' denotes 
the derivative 4>' = ^- F is an energy functional, whose hrst term Fi = dz gives the 
potential energy of the electric held over the insulator, and the second term F 2 = d^ 0^ 
the interaction energy between the electrostatic potential and charge stored in the device. 


A. Model device 

We consider a prototype memristor that contains donor defects with the activation energy 
Ed measured from the conduction band (CB). The insulator storage medium forms Schottky 
connections with both top and bottom electrodes with Schottky barriers ecpa or e06, and 
is assumed smaller than either of these barriers. Depending on initialization, we hnd three 
different families of numerical solutions to Eq. ([1]) illustrated by the band bending diagrams 
shown in Fig. 1 (a)-(c). The LRS shown in Fig. 1 (a) is characterized by a convex 0(z) 
where most of the CB is hlled by electrons and most donor defects are not ionized; because 
of the CB hlling, this state exhibits excellent band conductivity. The other convex 0(z) 
solution shown in Fig. 1 (b) is an IRS, which has the Fermi level {Ep) mostly pinned by 
Ep. Here most donor defects are not ionized, but because Ep is very close to Ep, we can 
expect quite large hopping conductivity. Note that if Ed is very close to the CB and there 
is signihcant overlap of Ep and CB at room temperature, the IRS becomes engulfed by the 
LRS, i.e., such IRS only appear for deep defects with Ep well inside the band gap. At the 
same time, in the case of multiple deep levels, there may be multiple IRS’s. Finally, the HRS 
shown in Fig. 1 (c) exhibits a concave 0(z) such that all donor defects are ionized, the CB 
is completely empty, and there are no carrier or impurity states nearby Ep. We may expect 
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FIG. 1. Band bending diagrams for different resistive states, (a) Low resistivity state LRS; (b) 
Intermediate resistivity state IRS; (c) High resistivity state HRS; (d) LRS at negative bias voltage 
V; (e) HRS at positive V. The strong bent lines indicate conduction and valence bands CB and 
VB and the dotted line indicates the donor activation energy Ed measured from CB. The dashed 
thin line indicates the Fermi level Ep. V applied to the top electrode shifts the band edges at 
z = a hy eV up for negative and down for positive V, as shown in (d) and (e). 


poor hopping conductivity at most, i.e., this state is an insulator. These different conductive 


mechanisms are supported by recent experimental observations of band conductivity for LRS 
and hopping conductivity for HRS [30, l3l|, ensued by large on/off resistivity ratios of 10^ 


or more. 


B. Resistive switching 

Applying a bias voltage of V to the top electrode shifts the band edges ai z = ahj —eV. 
Thus applying a negative bias in the LRS reduces the number of electrons in the CB and 
ionizes more and more donor defects, as shown in Fig. 1 (d), until the system exhibits a 
transition to the HRS. The transition from IRS to HRS occurs similarly. Applying a positive 
bias to the HRS, on the other hand, when V > (pa — EdIg, causes an increasing number of 
donor levels to be de-ionized, causing a transition to the IRS, and once V > (pa the CB 
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(b) 



Bias voltage V (V) 



Bias voltage V (V) 


FIG. 2. Energy F and charge Q stored in prototype memristor as a function of bias voltage V. 
(a) and (b) show the energy F and charge Q for defect scenario (I), a homogeneous distribution of 
oxygen vacancies Vq with concentration 0.005 nm“^. (c) and (d) show F and Q for defect scenario 
(II), a homogeneous distribution of oxygen vacancies Vo with concentration 0.005 nm“^ and Ti 
interstitals Th with concentration 0.0025nm“^. F and Q corresponding to HRS, IRS, and LRS are 
given by red, orange, and blue lines, and discontinuous jumps at switching between these states 
is given by the black line segments. Switching between HRS and LRS is indicated by the arrows 
labeled by 1 and 2, and switching between HRS and IRS by the arrows labeled 3 and 4. The insert 
in (d) shows a magnification of the HRS and IRS switching. 


becomes occupied, causing a transition to the LRS. It is important to note that while the 
above discussion holds for a system with mainly donor defects, the equations are symmetric 
with respect to charge, and, the above hndings hold equally well for systems with acceptor 
defects and hole carriers. 

To illustrate the switching mechanism, we calculate the energy F and charge Q of our 
prototype memristor as a function of bias voltage V at room temperature (T = 300 K), 
shown in Fig. 2. The red, orange, and blue lines indicate the HRS, IRS, and LRS, and the 
switching from one state to the other is indicated by arrows. Our prototype memristor is 
composed of a 20 nm slab of rutile Ti 02 sandwiched between metal electrodes that form a 
Schottky barrier of e0a = e0b = 1 eV at both sides. Here, rutile has a band gap of 3 eV, 
a relative permittivity of e = 80, effective electron masses of m* = 8.5 and 0.4 for heavy 
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and light electrons, and the effective hole mass of m* = 0.9 (effective masses were calcnlated 


from hrst principles and are well in agreement with experiment 


32|). We consider two defect 


scenarios that mimic the oxygen dehciency of Ti02 memristors: (I) rntile contains oxygen 
vacancies Vq with a concentration of 0.005 nm“^ (= 5 x 10^^ cm“^), which act as shallow 


donble donors 


33| : and (II) rntile contains 0.005 nm ^ of Vq together with 0.0025nm ^ of 


Ti interstitals Tij, which act as a qnac 


Ed = 0.6 eV below CB (see Refs. 


rnple deep donor with the donor activation energy at 
35| : also the Magneli phase hlament like structnres 


34 |35| 

rfl). 


observed in Ti02 may act as donors [36|). The memristor energy F for scenario (I) shown 
in Fig. 2 (a) has two energy branches, and for scenario (II) shown in Fig. 2 (c) has three 
energy branches, corresponding to the different resistivity states. These energy branches 
coexist over a wide range of V, and, most importantly, coexist at R = 0. The npshot is that 
once the system is in one of the different resisitivity states, it remains in that state nnless a 
bias in excess of the switching thresholds is applied. 


C. Charge and electric field inside memristor 


Figs. 2 (b) and (d) show the charge Q stored in the memristor. Notice that none of 
the states HRS, IRS, or LRS is charge nentral, i.e., Q ^ 0 always. The insnlator storage 
medinm of a memristor thns is spontaneonsly charged, snch that the HRS and LRS have 
opposite polarities. In onr example cases, for R = 0, 0 is symmetric, bnt in experimental 
devices, the Schottky barriers of the top and bottom electrodes are not necessarily eqnal, 
and the defect concentration is likely inhomogeneons. Snch asymmetries make the system 


polar, which combined with charge Q ^ 0 stored in t 


le device means that it can also act 


as a nanobattery, as observed in recent experiment When the HRS or IRS is switched 


to the LRS, there is a flow of electrons into the system, becanse the connection with the 
electrode is essentially ohmic (see Fig.l (e)). As R then decreases below the switching 
threshold, a Schottky barrier is formed again, preventing leakage of Q from storage medinm 
to electrodes. However, as R is farther decreased, Q gradnally approaches 0, and once the 
system is depleted from electrons, it is reset to the HRS; for the LRS, Q < 0, and for the 
HRS, Q > 0. The switching between HRS and IRS is similar, albeit here for both states 
Q > 0. As R —>■ (l)a — E^/e, electrons may freely occnpy the deep levels E^ dne to Tij, 
which canses a discontinnons drop in Q magnified in the inset of Fig. 2 (d). Switching back 
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FIG. 3. Electrostatic potential </> and electric field E inside memristor storage medium as a function 
of z for the 20 nm wide prototype memristor with defect scenario (II) at different bias voltages V. 
Notice that 4> is the electrostatic potential for electrons, i.e., the sign convention is opposite to the 
usual; with this sign convention (p illustrates the shape of the storage medium conduction band, 
(a) F = 0; (b) E = —0.3 V; (c) V = —1.9 V; (d) V = 0.4 V; and (e) V = 0.9 V. Red, orange, and 
blue lines correspond to HRS, IRS, and LRS. The horizontal dashed line indicates the constant 
electric held Ec = VjL. 


to HRS occurs as all electrons are depleted from the impurity states. 

We further illustrate the charging of the memristor by analyzing the electric held E 
inside the memristor storage medium. Fig. 3 shows the electrostatic potential (j){z) and 
corresponding electric held E = dcp/dz (there is no — sign because 0 denotes the electrostatic 
potential for negative charges; see App. for our prototype memristor in defect scenario 
(II) at zero bias [Fig. 3 (a)] and diherent bias voltages V immediately below switching 
threshold values. At zero bias for the HRS, > 0 at the top electrode and < 0 at the 
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bottom electrode, and thus any electrons possibly excited to the CB are expelled to either 
electrode. This combined with the fact that all donor defects are ionized gives Q > 0. For 
the LRS and IRS, electric held direction has changed and < 0 at the top electrode and 
£' > 0 at the bottom electrode, i.e., electrons excited to the CB will be trapped in the central 
region of the storage medium. Notice that for the LRS and IRS, most donor defects are not 
ionized, but nonetheless there are some electrons in the CB of the IRS, and for the LRS the 
CB is mostly hlled by electrons, and therefore Q < 0. Here it is important to recall that at 
both electrodes, there is a metallic reservoir of electrons, so any electrons stored inside the 
memristor (for LRS and IRS) are taken from these reservoirs, and vice versa for the HRS. 

The electric helds inside the memristor storage medium are very large—up to ±100 MV/m- 
even for zero bias voltage [Fig. 3 (a)]. These internal electric helds are of the same order of 
magnitude as constant electric helds Ec = VjL typically assumed to cause the accelerated 


dihusion 2J, [251, l27j of vacancies or impurity atoms and the ensuing formation or breaking 


of conductive hlaments. Comparison oi = V/L for large bias voltages around switching 
thresholds [dashed lines in Figs 3 (b)-(e)] leads to a dilemma: if an electric held of E >50 
MV/m would cause oxygen ions from the insulator storage medium to leave their lattice 
sites (form vacancies) and drift out of the material, or metal atoms from the electrodes to 
drift inside the storage medium, such a process would always take place for the convex band 
bendings (LRS or IRS) close to the electrodes, where the intrinsic electric held is larger than 
this even at zero bias. In fact, the band bending and intrinsic electric held at these elec¬ 


trodes is similar to what one obtains for a Schottky connection of the same materials 


, so 


applying the theory of accelerated dihusion [2J, l25|, l27| to a simple metal-insulator connec¬ 
tion would predict spontaneous ion dihusion and structural rearrangement within the space 
charge region. Thus, in order to understand the possible ion drift mechanisms or hlament 
formation mechanisms in memristors or other metal-insulator systems, one cannot rely on 
such simplihed approximations of constant electric helds {Ec ~ A(j)/Az)] instead, electric 
helds inside dielectric medium must be correctly evaluated as E = (—)V0 (see above for 
sign convention). 

Previous electronic switching theories suher from the voltage-time dilemma j^, but this 
is not the case for the present theory. This is easily illustrated by either Q or E inside the 
memristor storage medium. We must acknowledge that, in the LRS state, there is always a 
hnite probability Pout(i) for one electron to spontaneously tunnel out of the storage medium. 
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which for previous electronic models would cause the system to spontaneously forget its 
resistivity. In the present theory, however, LRS is a stationary state, and the probability 
for the next electron to tunnel out Pout( 2 ) is smaller, i.e., Pout( 2 ) < Pout(i)- By the same 
token, the probability for an electron to tunnel back into the storage medium from one of 
the electrodes Pin is larger than Pout( 2 ), i-e., the metastable state with one electron removed 
from the LRS tends to relax back into LRS rather than to decay into HRS. The same can 
be concluded based on the internal electric field in the storage medium: Figs 3 (b) and (d) 
show the electric fields for small bias voltages below switching thresholds. For LRS and IRS 
P < 0 at the top electrode and P > 0 at the bottom electrode, which maintains Q < 0, and 
vice versa for the HRS. After such ’’read” operations at zero bias, the intrinsic electric fields 
drive the system towards the stationary states shown in Fig. 3 (a). The charge difference 
between HRS and IRS in our prototype, however, is smaller than 1 electron per nm^, so the 
occasional tunneling of just one electron out of the system can cause the system to forget 
its state. Thus, for a reliable memristor, the difference in Q between the different resistivity 
states should be as large as possible. This charge difference can be maximized by choosing 
storage medium materials with a large electron (or hole) effective mass. 


III. MEMRISTOR MATERIALS 


The properties of the different resistivity states depend on the memristor material pa¬ 
rameters, defect abundances, and storage medium thickness L. An obvious prerequisite for 
resistive switching is the coexistence of the different resistive states. We investigate this 
coexistence by evaluating the energy functional P and its components Pi and P 2 as a func¬ 
tion of L for HRS, IRS, and LRS in two completely different materials, shown in Fig. 4. 
P and P 2 are shown for rutile Ti02 with Vq at concentrations of 0.001 nm“^ [Fig. 4 (a)], 
0.005 nm“^ [Fig. 4 (b)], and 0.01 nm“^ [Fig. 4 (c)] and a correlated concentration of Ti, 
([Tii] = |[Vo]). The same are shown for ZnO with Cr impurities at the concentrations of 
0.0001 nm“^ [Fig. 4 (d)], 0.0005 nm“^ [Fig. 4 (e)], and 0.001 nm“^ [Fig. 4 (f)]. Cr in ZnO 
is a deep impurity with two donor levels at 1.0 and 2.6 eV below the conduction band [1 


ZnO has a band gap of 3.4 eV, effective masses m* = 0.28 and m* = 0.6, and a relative 
permittivity of ~ 10. P for the HRS increases rapidly as L increases, mainly due to the 
similar increase of P 2 . At P larger than some threshold value, P for HRS becomes exceed- 
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Storage media thickness L (nm) 


FIG. 4. Memristor energies F as a function of storage medium thickness L. (a), (b), and (c) show 
F and its component F 2 for rutile Ti02 with Vo at the concentrations of (a) 0.001 nm“^, (b) 
0.005 nm“^, and (c) 0.01 nm“^ and a correlated concentration of Tii ([Tii] = ^[Vq])- (d), (e), 
and (f) show F and F 2 for ZnO with Cr impurities at the concentrations of (d) 0.0001 nm“^, (e) 
0.0005 nm“^, and (f) 0.001 nm“^. The red, orange, and blue solid lines show HRS, IRS, and LRS, 
and the dashed lines show the corresponding F 2 . 


ingly large >> j ^ ^ HRS becomes unstable. This sets a maximum size for the 

memristor, which for Ti02 and ZnO is estimated to be some tens of nm depending on the 
defect concentration. 

The superlinear behavior of F 2 described above arises as follows. F 2 is proportional to 
p multiplied by 0 integrated over the length L of the memristor. For the HRS, all defects 
are ionized so p is constant, proportional to defect concentration Nd. At the same time, 
for HRS, 0 is increasingly bent upwards as L increases, which leads to the superlinear 
behavior of F 2 . The curvature of 0 is inversely proportional to the relative permittivity, and 
hence, F 2 increases much more rapidly in ZnO than Ti02 at similar defect concentrations 
(cf. panels (a) and (f) in Fig. 4). For small L (around 10 nm or smaller), HRS, IRS, 
and LRS are indistinguishable, because there is no band bending. This is because the 
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curvature of 0 is proportional to defect concentration. It follows that in ZnO, where F for 
the HRS is reasonably small only for very small defect concentrations, the IRS and LRS 
with convex band bendings only emerge at rather large L. The above discussion illustrates 
the interplay of relative permittivities and defect concentrations, and serves as a guideline 
to design memristors using arbitrary insulator materials as storage medium. 


A. The role of filamentation in resistive switching 


Filament formation inside the memristor storage medium has been observed in various 
devices 


14-18 


20|-|22| . For example, conical regions of oxygen dehcient Magneli phases have 


been observed in Ti 02 memristors, and chain like structures of metal impurity atoms have 


been observed in Conductive Bridging Random Access Memory (CBRAM) devices 401. l41| 


It is clear that many memristor devices contain such hlament like structures, but their role 
in resistive switching is not so obvious. While voltage controlled hlament formation and 
manipulation has been demonstrated, the time scales for such concerted ionic motion are 
much slower (seconds or minutes) than typical resistive switching times (nanoseconds or 
less). A recent in-situ transmission electron microscopy (TEM) study of a memristor device 
shows that resistive switching is accompanied by the appearance and disappearance of a 
dark region in TEM images ensuing the switching to LRS and to HRS, respectively, which 
was interpreted as the formation and erasure of conductive hlaments 23(]. However, these 
changes appear in the TEM images after the switching, over a time period of ~ 0.15 seconds, 
which would suggest that the contrast change in TEM is a consequence of the switching, 
rather than its cause. 

One of the present authors has shown that quasi-ID metal impurity chains inside a 
semiconductor medium behave akin to a large concentration of deep impurities the semicon¬ 
ductor rather than a metallic inclusion. These ID chains exhibit multiple (deep) donor 


or acceptor transitions, shifted from similar transitions observed for isolated impurities. We 


have also shown that the 
an intermediate band 36 


iO Magneli phases exhibit a semiconductor band structure with 


43|, which essentially may behave like a deep donor state. Thus, 


instead of the hlaments behaving like a metallic lead connecting the electrodes, we propose 
that their role is electronic, i.e., they behave as donor or acceptor type defects, i.e., in the 
LRS, the conductivity is not through the electronic states of the hlament, but the insulator 
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Magnetic configurations of LRS 
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FIG. 5. Spin configurations of magnetic memristor. The polarity of the bottom electrode is 
assumed to be fixed (red arrow), whereas the top electrode polarity (blue arrow) is assumed to be 
easily reversible (writable). The storage medium polarity (purple arrow) can be controlled by the 
polarity of the top electrode at the time of switching. 

conduction band. 


B. Magnetic memristor 


The above case study of a memristor based on Cr-doped ZnO discussed above leads to 
an intriguing extra feature. Transition metal doped insulators may exhibit carrier-mediated 
ferromagnetism 4^, l4^. For Cr-doped ZnO (or In 203 ) dy, the mechanism of magnetism 
is such that the ferromagnetic interaction is switched on by degenerate electron doping, 
i.e., by hlling the CB by carrier electrons. Thus, the LRS of a memristor using Cr-doped 
ZnO or In 203 would be spontaneously magnetized. The polarity of the storage medium 
magnetization could be controlled by a ferromagnetic top electrode, i.e., the storage medium 
is assumed to adopt the polarity of the top electrode at the instant of switching. Such a 
magnetic memristor has four spin conhgurations in the LRS, as shown in Fig. 5. For an 
asymmetric memristor, all four spin conhgurations are expected to have different resistivities 
due to spin-dependent transport (see e.g. Ref. [48|), leading to increased storage capacity 


of the memristor device. T 


half-metal 

interface. 


ic ferromagnet 


le spin-dependent transport could be enhanced by having a 


49 


m 


50| as one of the metal contacts, such as e.g. the GdN/GaN 
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IV. CONCLUSION 


We present an electronic mechanism of resistive switching for nanodevices. This mecha¬ 
nism is demonstrated by self-consistent calculations of the electrostatic potential (j){z) from 
Poisson’s equation for prototype resistive memory devices, composed of an insulator slab 
sandwiched between metallic electrodes that form Schottky barriers with the insulator. The 
two key ingredients for resistive switching are (i) the insulator must be abundant with donor 
or acceptor type defects and (ii) the insulator slab thickness has to be less than some tens 
or hundreds of nanometers; the maximal thickness depends on the insulator properties and 
defect abundance. When both conditions are met, the system exhibits multiple band bend¬ 
ing solutions, which may correspond to the high resistivity state, where the insulator is 
completely depleted of charge carriers, the low resistivity state where the insulator valence or 
conduction band contains a maximal amount of charge carriers, or intermediate resistivity 
states, where the Fermi level is pinned around deep impurity levels inside the gap. The right 
amount of defects can be obtained during initial electroforming, which has been demon¬ 


strated to form various kinds of defects in the insulator storage medium 


, 18 |, 


10 


18 


20 |. 


These different 0 are stationary states, robust against small perturbations, and stochastic 
tunneling of electrons into or out of the insulator storage medium cannot change the resis¬ 
tive state of the device—switching occurs only once the bias voltage exceeds given threshold 
values. Due to the electronic nature of the switching, it may be possible to switch resistive 
states also optically, provided the device is sufficiently optically transparent. This electronic 
nature also explains the fast switching times. Finally, we propose a magnetic memristor, 
which could allow simultaneous control of both 0 and magnetic states of the system, and 
increase the storage capacity of the memristor. 


Appendix A: Self-consistent calculations 

Let us consider an insulator material sandwiched between two metal electrodes at z = a 
and z = b, which form Schottky contacts with n type Schottky barriers e(j)a and e0b, i.e., a 
one-dimensional Poisson equation with the Dirichlet boundaries 0(a) = 0^ and 0(6) = 0;,. 
0(z) is the electrostatic potential of an electron (charge —e) chosen such that e0 marks the 
position of the conduction band minimum with respect to electron chemical potential (Fermi 
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level) fixed at e0 = 0. To find the 0 that minimize Eq. (| 

d^0 p (p dp 

dz"^ eeo eeo d(p 

self-consistently with the charge density p. 

The charge density in Eq. flAip is 


p = -He + Pv + A^D(gD,o - ^d) • 


we solve the Euler equation 

(Al) 


(A2) 


p is composed of the densities of electrons in the conduction band (ric), holes in the valence 
band (p^), and the concentration of donor or acceptor defects (iVo) multiplied by the charge 
of that defect gD,o when all its gap levels are empty minus the number of electrons occupying 
the defect induced levels (ud). Electron and hole densities in conduction and valence bands, 
He and py are calculated as 

9ci£) 


POO 

nc= d£ / 

Je(i> exp{e/kT) + 1 

p.= de 


(A3) 

(A4) 


exp(-£/A;T) + 1 ‘ 

£ is the electron quasiparticle energy, and and pc are effective densities of states of the 
valence and conduction bands given by 


t.3/2 


9ci.£) = ^|^ 72 (£-e 0 ) 


= 3/2 


Tfi I 

9v{^) — W3 2(e0 — £gap ~ ^) • 

n TT^ ’ 


(AS) 

(A6) 


Here m* and m* denote the effective masses of the conduction and valence bands, and h is 
the Planck constant. The number of electrons occupying donor or acceptor levels in the gap 
is given by 


riD = 


J2Njexp{-ej/kT) 


(A7) 


Eexp(-£j//i;T) 

where ej and Nj are the energy and number of electrons in state j, and the summations are 
taken over all electronic states within the gap. 

A self-consistent solution is obtained as follows. We start by evaluating p (Eq. (IA2p i for 
some trial (p. Then we solve Eq. flAip numerically, using the hnite difference approximation, 
to obtain a new (p. Using Eq. flA2p again, we construct a new p. We continue such iterations 
until the energy given by Eq. ([2]) converges. Starting from different trial (p can lead to 
completely different solutions. 
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SUPPROTING INFORMATION 

The resistive switching is demonstrated by a movie illustrating the evolution of the elec¬ 
trostatic potential cj){z) during voltage sweeps. 

TOC ENTRY 

Sandwich structure nanodevice has different band bendings that correspond 
to different resistivities. For high resistivity state, the insulator storage medium is 
completely depleted of charge carriers; for low resistivity state, storage medium contains 
maximal amount of charge carriers. Bias voltage shifts top electrode Schottky barrier, and 
at given threshold bias, system switches resistivity. 
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